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ABSTRACT 

The spectrum of the supernova relic neutrino background (SRN) from collapse-driven 
supernovae ever occurred in the universe is calculated by using a realistic, time-dependent 
supernova rate derived from a standard model of galaxy evolution based on the population 
synthesis method. The SRN spectrum we show here is the most realistic at present, because the 
largest uncertainty in previous theoretical predictions has come from unrealistic assumptions 
of the supernova rate so far made. The SRN is one of the targets of the Superkamiokande (SK) 
detector which will be constructed in a year and the SRN, if at all detected, would provide a 
new tool to probe the history of supernova explosions in the universe. The expected event rate 
at the SK is therefore calculated in this paper. Our major results include: (1) the supernova 
rate is much higher in the early phase of evolution of galaxies and there appears a hump in 
the SRN spectrum in the low-energy region of < 5 MeV, (2) the SRN flux depends on the 
Hubble constant (H ) in a way approximately proportional to Hq and only weakly on the 
density parameter of the universe (Slo) and a cosmological constant (Ao), (3) the uncertainty 
in the star formation history of spiral galaxies affects the resulting SRN flux by about a factor 
of 3, and (4) the plausible event rate at the SK is 1.2 yr _1 in the observable energy range 
of 15-40 MeV. Such a low event rate is due mainly to a quite low supernova rate at present 
which is averaged over the morphological types of galaxies. The most optimistic rate in our 
model is found to be 4.7 yr _1 in the same energy range, and if more events are detected, we 
will have to reconsider our current understanding of collapse-driven supernovae and evolution 
of galaxies. 

Subject headings: cosmology: diffuse radiation — galaxies: evolution — stars: supernovae: 
general 
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1. Introduction 

The water Cerenkov neutrino detector, Superkamiokande (SK) (Nakamura et al. 1994), will be 
constructed and observations will start in a year. One of the targets of this detector is to detect the 
supernova relic neutrino background (SRN) which would originate from all collapse-driven supernovae ever 
occurred in the universe (Bisnovatyi-Kogan k Seidov 1984; Krauss, Glashow, k Schramm 1984; Dar 1985; 
Woosley, Wilson, k Mayle 1986; Zhang et al. 1988; Hirata 1991; Totsuka 1992; Suzuki 1994; Totani k 
Sato 1995). Since the emergence of neutrinos from a supernova has been confirmed from both numerical 
simulations and observations of SN1987A, comparison of the SRN spectrum with observation at the SK will 
enable us to constrain a supernova rate and therefore provide us valuable information on the star formation 
history in galaxies. In this regard, in order to foresee a capability of the detector it is worth predicting 
the SRN spectrum and the resulting event rate at the SK for a given rate of supernovae in various models 
of the universe, either with or without a cosmological constant. In order to detect the weak SRN signal, 
various other backgrounds must be removed. Fortunately, there is a quite silent window between 15 and 40 
MeV where various backgrounds other than the SRN are strongly suppressed (Hirata 1991, Totsuka 1992). 
Although it is not certain whether the SRN can be detected even in this window, it is important to predict 

the most believable flux of the SRN which has only roughly been estimated in a broad range of 1-10 2 

— ? — 1 
cm s . 

A major source of the uncertainty in previous theoretical predictions has been a supernova rate 
prescribed in the models. While most of previous papers assume a supernova rate which is constant with 
time, it is more likely that the past rate is higher than at present as galaxies evolve. The effect of galaxy 
evolution is considered by Bisnovatyi-Kogan k Seidov (1984), assuming that a supernova rate depends 
on time or redshift in a manner proportional to (1 + z) 7 , where z is the redshift parameter. This effect 
is also considered in our previous paper (Totani k Sato 1995, henceforth Paper I), assuming that half of 
the supernovae have exploded like a burst at a certain redshift and another half have ever occurred at a 
constant rate. These assumptions are, however, rather artificial and a reasonable model of galaxy evolution 
should be used for prescribing a realistic, time-dependent supernova rate. 

In this paper, using a standard model of evolution of galaxies based on the method of stellar population 
synthesis (Arimoto k Yoshii 1986, 1987; Arimoto, Yoshii, k Takahara 1992), we calculate a supernova rate 
as a function of look-back time of galaxies and predict the most realistic SRN spectrum as well as the 
event rate at the SK for various cosmological models. This model of galaxy evolution has been used to 
investigate the geometry of the universe from comparison with the number counts of faint galaxies (e.g., 
Yoshii k Takahara 1988; Fukugita et al. 1990; Yoshii 1993; Yoshii k Peterson 1995) and basically agrees 
with observations, particularly after the selection effects are taken into account in the theoretical framework 
(Yoshii 1993). Some variants of the standard model of galaxy evolution are also used to delineate a range 
of uncertainty in the resulting SRN spectrum arising from different histories of star formation in galaxies. 

Throughout this paper, we consider only electron antineutrinos (jVs) from collapse-driven supernovae, 
because the contribution of neutrinos from thermonuclear supernovae to the SRN is negligible and i> e 's are 
most easily detected in a water Cerenkov detector. In §2, we construct a model of evolution of galaxies, the 
supernova rate, and neutrino emission from a supernova. Formulations for calculating the SRN spectrum 
are given in §3. The results are given in §4, and their dependence on models of the universe as well as 
galaxy evolution is also discussed. After the total (energy-integrated) SRN flux of our results is compared 
with previous estimates in this section, we make a prediction for the response of the SK detector in §5, and 
finally summarize the results in §6. 
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2. Model Construction 



2.1. The Population Synthesis Model of Galaxy Evolution 

We summarize here the model of evolution of galaxies used in this paper, following Arimoto & Yoshii 
(1986, 1987) and Arimoto, Yoshii, & Takahara (1992). It is assumed that a galaxy consists of three 
components, namely, interstellar gas, stars, and stellar remnants. Time evolution of the mass fraction of 
gas in a galaxy, f(t), is determined by the differential equation: 

f t =-C(t) + F(t) + A(t), (1) 

where C(t) is the star formation rate (SFR), F(t) is the ejection rate of gas from dying stars, A(t) is the gas 
accretion rate (ACR) from the surrounding region of a galaxy, and t is the time elapsed from the formation 
of a galaxy, not a cosmic time. 

It is generally assumed that the star formation law can be approximated as 

C(t) = „/»(*), (2) 

where v is the SFR coefficient and has the dimension of (time) -1 , and n is the power index. A large 
number of papers constraining the form of SFR from observations show that n is most likely between 1 and 
2 for spiral galaxies (e.g., Kennicutt 1989), and we consider two cases of n = 1 and n = 2 in this paper. 
For elliptical galaxies, however, the remaining gas would leave the system by a galactic wind driven by a 
burst of supernova explosions in the early stage of evolution of galaxies. We introduce this effect in the star 
formation law for elliptical galaxies as follows (Arimoto & Yoshii 1987): 

c(t)={: m (3) 

L U (t > t GW ) , 

where taw is the time when the galactic wind occurs. 
The ejection rate, F(t), is expressed as 



F(t)= C{t -tm^mRmdm . (4) 

J mi 

Here, m is the stellar mass, R m is the fraction of ejected mass from a dying star, <f> m is the initial stellar 
mass function (IMF), and m u and m; are the upper and lower bounds of the stellar mass, respectively. In 
this paper we adopt mj = O.O5M and ra„ = 6OM . The lifetime of a star of mass m is denoted by t m and 
is taken as t m = 11.7(m/M0) _2 Gyr (Talbot & Arnett 1971). The IMF <f> m is assumed to be of a power-law 
form: 



(ft - l)mf 1 



(5) 



1 - (rrn/m u y 1 

where \i is the power index and the function is normalized to unity. The fraction of ejected mass from a 
dying star is assumed to be 

R _ ( (m, <m < m rem ) 

m \ l-(m rem /m) (m rem < m < m u ) , 

with m rem set to be O.7M (Talbot & Arnett 1971; Arnett 1978). 
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Solutions to the above equations depend on the assumptions made for gas accretion from the 
surrounding region of a galaxy. One is the 'simple', one-zone model, in which there is no accreting gas 
from the surrounding region, and all the gas contributes to star formation at the initial stage of galaxy 
evolution, i.e., A(t) = and /(0) = 1. The other is the 'infall' model, in which there is no gas initially 
and the gas is supplied gradually from the surrounding region of a galaxy. The form of ACR is assumed 
as A(t) = aexp(— at) and /(0) = (Arimoto et al. 1992), where a is the ACR coefficient and has the 
dimension of (time) -1 . Chemical evolution studies suggest that gas infall plays a decisive role in reproducing 
the photometric and chemical properties of present-day spiral galaxies (e.g., Matteucci & Francois 1989; 
Arimoto & Jablonka 1991). We therefore consider two alternative assumptions such that the accretion is 
turned on and off for spiral galaxies. On the other hand, for elliptical galaxies, we consider only the simple 
model without gas infall. 

In this paper, we group E and SO galaxies together and classify spirals into four subclasses like 
Sab, Sbc, Scd, and Sdm. Evolutionary behaviors of chemical abundances in the gas and stars as well as 
integrated colors of stars in a galaxy are modeled with the parameters of v, taw an d H for ellipticals, or 
those of v, a, n and \i for spirals. These parameters have been fixed through a fit to the observed colors of 
galaxies. For spirals from Sab to Sdm, the adopted values off and a for n = 1 or 2, together with Salpeter's 
IMF (jjl = 1.35), are the same as those in Arimoto et al. (1992). For ellipticals, the adopted values of v 
and a corresponding taw are those for the lower, standard, and higher SFR models in Arimoto & Yoshii 
(1987); the lower and higher bounds of v correspond to an intrinsic dispersion of SFR, which can account 
for the observed scatter along the mean relationship on the color-magnitude (CM) diagram. We note that a 
somewhat flat IMF with \i = 0.95 well reproduces the observed optical through infrared colors of ellipticals 
(Arimoto & Yoshii 1987) and we adopt this value of \i for E/SO-type galaxies. 

In this paper, we consider four cases (simple /infall and n = 1 or 2) for each of Sab, Sbc, Scd and Sdm, 
and three cases (lower, standard, and higher SFRs) for E/SO-type galaxies; these models are considered to 
cover the uncertainty of the evolution model of galaxies. Except when the uncertainty of galaxy evolution is 
concerned, we use the SI model to represent the evolution of spiral galaxies, and the standard SFR model 
for E/SO-type galaxies; here and hereafter, for simplicity, we refer SI or S2 to the simple model with n = 1 
or 2, and similarly II or 12 to the infall model with n = 1 or 2. 



2.2. Comparison with the Number Counts of Faint Galaxies 

In a number of papers trying to determine the geometry of the universe from the number counts 
of faint galaxies, the above model of galaxy evolution has been used and basically shown to agree with 
observations in some different wave bands (e.g., Yoshii & Takahara 1988; Fukugita et al. 1990; Yoshii 1993; 
Yoshii & Peterson 1995), while the values of the density parameter of the universe (Slo) and a cosmological 
constant (A ) have not yet been settled. In this section, we compare the predicted number counts of faint 
galaxies with the observations to show the validity of using the model constructed in the previous section. 

The if -band data of number counts are those compiled in Yoshii & Peterson (1995) with the addition 
of some new observations (Djorgovski et al. 1995; McLeod et al. 1995; Aragon-Salamanca & Ellis 1994). 
Among all these surveys, Djorgovski et al.'s survey uses the faintest surface brightness threshold for image 
detection, thus providing the faintest number count data. Given a fixed threshold, the selection effects 
become quite significant near the survey limit (Yoshii 1993) and these effects should be taken into account in 
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the predictions to be compared with the uncorrected counts in the faint end. The observational conditions 
used to predict the uncorrected counts are matched with those actually used by Djorgovski et al. (1995), 
that is, the 1.5" aperture magnitude with Sl = 24 mag/arcsec 2 and the typical seeing of 1" FWHM, where 
Sl is the threshold of surface brightness. For the minimum angular diameter, we adopt D m i„ = 1". The 
manner of including the selection effects is described in Yoshii (1993). 

Figure 1 shows the observations of the A'-band number counts and the predicted curves from the 
standard model of galaxy evolution (the Sl model is used for evolution of spiral galaxies, and the standard 
SFR model for ellipticals). The data are taken from various surveys with different SVs and the error bars 
represent the Poissonian error. The data with brighter Sl's should be corrected for undetected images 
when plotted in the same figure with Djorgovski et al.'s data with the faintest Sl- Therefore, except for 
Djorgovski et al.'s uncorrected counts to which the predictions are compared, shown are the corrected 
counts from all the other surveys with brighter SVs. 

Three different cosmological models are considered: a low density, flat universe with Slo = 0.2, Ao = 0.8, 
and h = 0.8 (thick line), a low density, open universe with Sl = 0.2, A = 0.0, and h = 0.8 (dashed 
line), and an Einstein-de Sitter universe with Slo = 1.0, Ao = 0.0, and h = 0.5 (dot-dashed line), where 
h = i?o/(100km/s/Mpc). Note that changing the Hubble constant hardly affects the theoretical predictions 
of number counts. The parameters of the luminosity function are those in Efstathiou, Ellis & Peterson 
(1988). This figure shows that the model predictions are well consistent with observations, while the model 
of the universe cannot be determined strictly. We calculate the SRN spectrum with these three cosmological 
models. 

Figure 2 shows the observations of the number counts and theoretical predictions of the model of 
galaxy evolution for only a low density, flat universe (Slo = 0.2, Ao = 0.8, h = 0.8) with the four different 
models of evolution of spirals; Sl (thick line), S2 (short dashed line), II (dot and short-dashed line), and 12 
(dot and long-dashed line). The effects of luminosity evolution are least in redder passbands because the 
steady accumulation of long-lived, low-mass stars in galaxies makes the evolution of the red light secular 
and weak (Bruzual 1983; Arimoto & Yoshii 1986), and all of these four models are well consistent with the 
observations of the number counts of faint galaxies. 



We have constructed the model of evolution of galaxies from which the supernova rate is obtained. 
The neutrinos of our concern are emitted from collapse-driven supernovae having massive progenitor stars, 
and their rate per unit mass of a galaxy is expressed as 



where m t is the mass for which t m = t, and msNi is the critical mass above which all stars end their life 
with collapse-driven supernova explosions; we adopt msNi = 8M throughout this paper. The SFR per 
unit mass of a galaxy, C(t), is numerically calculated from equations in §2.1. 

The supernova rates thus calculated are shown in Figs. 3 to 5, in units of yr _1 in a galaxy with the 
typical mass of IO^Mq. Figure 3 shows by the dashed lines the supernova rates for five types of galaxies 
using the Sl model for spirals and the standard SFR model for ellipticals. The supernova rate averaged over 



2.3. Supernova Rate 




(0 <t <t m J 
(j> m m~ 1 C(t -t m )dm (t>t mu ) , 



(7) 
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the types with a weight of a^(M/ Ls)k is shown by the thick line, where is the relative proportion of the 
k-th type galaxies (Pence 1976) and (M/Lb)ic is their 5-band mass-to-luminosity ratio (Faber & Gallagher 
1979). The values of and (M / ' Lb)^ are tabulated in Table 1. We note that the 'mass' of a galaxy in this 
paper represents the mass of the ordinary matter, i.e., field stars, remnants, and interstellar gas which is 
used in star formation. The mass-to-luminosity ratio is that within the Holmberg radius, within which the 
luminous and ordinary matter is considered to be dominant. Figure 4 shows the weighted average rates 
with the four different models of SI, S2, II, and 12 for spirals provided that the standard SFR model for 
ellipticals is used in common. Figure 5 shows the weighted average rates with the three different models 
(lower, standard, and higher SFR) for ellipticals provided that the SI model for spirals is used in common. 

In all the cases, supernovae most dominantly occur in the early phase of galaxy evolution ( < 1 Gyr) 
in E/SO-type galaxies, and then in the later phase ( > 1 Gyr) in spirals. This is a general trend largely 
irrespective of different assumptions made for galaxy evolution, as clearly seen in Figs. 3 to 5. It should 
be noted that the weighted average rate in the early phase is about two orders of magnitude higher when 
compared to that in the recent past, and more than the half of a total number of supernovae have exploded 
during the initial 1 Gyr; this feature affects strongly the shape of the resulting SRN spectrum. 

Now we compare our predicted average rate of collapse-driven supernovae with observations. There 
are various observational estimates of the supernova rate, and they are scattered with a mean around one 
collapse-driven supernova in 50 years in our Galaxy (Tammann et al. 1994). In general, our Galaxy is 
considered to be an Sbc galaxy with an age of 10-15 Gyrs, having a total mass of about IO^Mq. With 
these typical values of age and mass, our predicted supernova rate for Sbc spans from one in 50 years to 
one in 100 years (Fig. 3) in fairly good agreement with the observational estimate. 

We also examine whether or not our model meets the nucleosynthesis requirement. Since the oxygen 
is synthesized in the interior of massive stars and intermediate mass nuclei burn away to create iron in 
thermonuclear supernovae, the amount of oxygen is related directly to the time-integrated number of 
collapse-driven supernovae. Figure 6 shows the time-integrated number of collapse-driven supernovae in 
a galaxy with the mass of 10 11 M e as a function of time from the formation of the galaxy; results for 
individual types of galaxies are shown by the dashed lines and those for the weighted average by the thick 
line. Here, as in Figure 3, the SI model is used for spirals and the standard SFR model for ellipticals, but 
the results are hardly changed by other choices of evolution models. Taking the mass fraction of oxygen in 
the Galactic disk as about 8.5 x 10 -3 from observations of galactic HII regions (e.g., Diaz 1989), the total 
mass of oxygen in a galaxy of 10 n M Q is estimated as 8.5 x 10 8 M Q . Arnett, Schramm, & Truran (1989) 
showed that about 10 9 collapse-driven supernovae should be responsible for the oxygen enrichment to that 
amount in our Galaxy. Our predicted number for an Sbc galaxy is about 10 9 at present, which agrees with 
that required from the nucleosynthesis arguments. This agreement makes the evolution models used in this 
paper more believable, because the models have been contrived to reproduce the observed colors of galaxies, 
but not the cumulative number of supernovae in galaxies. 



2.4. Electron Antineutrino Emission from a Supernova 

Emission of electron antineutrinos (jVs) from a supernova is treated in the same way as Paper I. We 
assume that the energy distribution of the neutrinos obeys the Fermi-Dirac distribution without chemical 
potential. Numerical calculations of Woosley et al. (1986) for gravitational collapses of massive stars with 
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10, 15, and 25M are used to represent the total energy and temperature of i> e 's from a supernova in the 
three intervals of stellar mass such as 8-12. 5Mq, 12.5-2OM , and 2O-6OM , respectively. Contrary to 
Paper I where the number of stars in these three intervals is assumed to be constant with time, the number 
of stars changes with time in a way consistent with the model of galaxy evolution considered in this paper. 
Following the above approximation, we denote a differential number of emitted neutrinos from a supernova 
as dN v (q) / dq where q is the neutrino energy. 



3. Formulations 



The differential number flux of SRN is expressed as dF v (q) / dq = cdn v (q) / dq where dn v (q) / dq is the 
present number density of SRN per unit energy of neutrinos, and c is the speed of light. The neutrinos 
emitted in the interval of cosmic time between t and t + dt contribute to dn v (q) as follows: 



dn v (q) = p G (t)RsN(t ~ t F )dt 



dN v {(l + z)q} 



{(1 + z)dq}(\ + z)- 



(8) 



where pa(t) is the mass density of galaxies in the universe, t F is the cosmic time when galaxies were formed, 
z is the redshift parameter, and Rsn(^ ~ t F ) is the supernova rate at the cosmic time t per unit mass of 
galaxies. Note that the definition of Rsn in this paper is different from that in Paper I. The differential 
number of emitted neutrinos from a supernovae, dN v (q) / 'dq, is given in §2.4. The factor (1 + z) -3 comes 



from the expansion of the universe, and using pa(t) 

dn v (q) 



/°g(^o)(1 + z) , equation (8) reduces to 



Pcj{to)RsN{t ~ t F ) 



dN v {(l + z)q} 



(1 + z)dt 



(9) 



In order to proceed further, the quantity PG(to)RsN(t ~ t F ) is given by 

' M 



PG(to)RsN(t — t F ) 



E 

k 



dLLipk(L) 



L 



Rk(t — tF) 



(10) 



where ip^^L), (M / ' L)^, and Rk(t — tp) are the luminosity function, the mass-to-luminosity ratio, and the 
supernova rate for the k-th type galaxies, respectively, and the subscript k runs over five types of galaxies 
(E/SO, Sab, Sbc, Scd, Sdm). By definition, tpk(L) is the number density of the k-th type galaxies having 
the absolute luminosity L. We assume that the shape of the luminosity function does not depend on the 
type, yielding tpk(L) = aktp(L) where is the relative proportion of each type and tp(L) is the luminosity 
function for the composite of all types of galaxies. Furthermore, we assume that the mass-to-luminosity 
ratio does not depend on the luminosity of galaxies. Accordingly, equation (9) becomes 



dn v (q) 



C 



Rk(t — t F ) 



dN v {(l + z)q} 



(1 + z)dt 



(11) 



where C = f dLLi\>{L) is the luminosity density of galaxies. Integrating over z and using a relation 
dF v (q) I dq = cdn v (q) / dq, we finally obtain the differential number flux of SRN: 



dF v (q) 



cC 



Rk(t — t F ) 



dN v {(l + z)q} 



(1 + *) 



dz 



with 



dz 



-H (l + z)y/{l + n z)(l + zf - A (2z + z 2 ) 



(12) 



(13) 
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where zp is the formation redshift of galaxies, and Ao is the normalized cosmological constant; 
Ao = (c 2 A)/(3i?o 2 )- The time t — tp is related to z via t — tp = f* F \dt/dz\dz. Throughout this 
paper, following Efstathiou, Ellis & Peterson (1988), we adopt the 5-band luminosity density of 
C B = 1.93to:e x 1Q8/l L bq Mpc" 3 , where h = ifoAlOOkms^Mpc" 1 ). The recent results of the luminosity 
function (Loveday et al. 1992; Marzke et al. 1994) are also consistent with this value of Cp- The adopted 
values of (Pence 1976) and {M/LB)k (Faber & Gallagher 1979) are summarized in Table 1. With the 
above values of a^, (M / L)^, and C, the contribution of the mass of galaxies to the density parameter of the 
universe, Q,fj, becomes about 0.7%. We set zp = 5, except for otherwise studying the effect of changing zp 
in §4.2. 

In most of previous papers including Paper I, the absolute number of supernovae in the model 
is determined by a fit to the observation of the present supernova rate in our Galaxy, which is quite 
uncertain. In this paper, instead of using the present supernova rate, however, we use the observations of 
the luminosity density of galaxies in the universe and the mass-to-luminosity ratio of galaxies, and then 
derive the absolute number of supernovae in the model. 



4. Supernova Relic Neutrino Background (SRN) 

4.1. SRN Spectrum in Some Cosmological Models 

We calculate the SRN spectrum in some models of the universe, using the SI model for evolution 
of spirals and the standard SFR model for ellipticals. We use a low-density, flat universe with Slo = 0.2, 
Ao = 0.8, and h = 0.8 as a standard cosmological model. The SRN spectrum for this standard cosmological 
model and two variants with h = 0.5 and 1.0 are shown in Figure 7, by the thick lines. The spectra for a 
low-density, open universe with Sl = 0.2, A = 0.0, and h = 0.8 (dashed line) and an Einstein-de Sitter 
universe with Slo = 1-0, Ao = 0.0, and h = 0.5 (dot-dashed line) are also shown in this figure. 

For a purpose of comparison, the SRN spectrum with a constant supernova rate (from Paper I) is also 
shown by the dotted line for a standard cosmological model (Slo = 0.2, Ao = 0.8, h = 0.8). In order not to 
violate the nucleosynthesis requirement, the value of the constant supernova rate (1.6 x lO'V^Mpc" 3 ) 
is chosen such that the total SRN flux J (dF v j 'dq)dq, or equivalently the total number of supernovae ever 
occurred, is the same as that for the time-dependent supernova rate considered in this paper. Thus, from 
these two SRN spectra we clearly see how the effect of galactic evolution modifies the resulting SRN 
spectrum. 

In all SRN spectra, there appears a hump in the low-energy region of < 5 MeV. This feature is similar 
even to the constant+burst model in Paper I. The reason is because a very large number of neutrinos 
emitted from supernova bursts in the early phase of galaxy evolution are strongly redshifted to populate 
the low-energy part of the spectrum. On the other hand, in the high-energy region of > 10 MeV, most of 
neutrinos are from supernovae which occurred at a smaller rate in the recent past. 

The SRN flux depends most significantly on the Hubble constant among the cosmological parameters, 
and it approximately scales as dF v / dq oc Hq . This is understood as follows: The SRN flux is proportional to 
the mass density of galaxies in the universe, M, and its observational estimates are generally proportional 
to Hq like M = C x (M/L) oc H x H . Furthermore, the total number of supernovae and therefore the 
total SRN flux apparently increase with increasing the age of galaxies, or equivalently with decreasing Hq 
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because the age is proportional to Hq . This effect is, however, very small; most of the supernovae have 
exploded within the initial 1 Gyr from the formation of galaxies, so that the total number of supernovae 
does not significantly change with the age of galaxies (Fig. 6). On the other hand, shortening of the age 
leads to a higher present supernova rate (Fig. 3), therefore the SRN flux is enhanced in the high-energy 
region ( > 10 MeV) where neutrinos are easily detected in the SK detector. As a result, the event rate at 
the SK scales as Hq or a power of Ho with a little higher index. 

Within a considered range of Slo = 0.2-1, Ao = 0-1, and h = 0.5-1, the SRN spectrum is less sensitive 
to Slo or Ao when compared to Ho- Smaller Slo or larger Ao decreases the SRN flux in the high-energy 
region, since these effects equally prolong the age of galaxies and therefore decrease the supernova rate in 
the recent past. 



4.2. Uncertainty from Models of Galaxy Evolution 

We estimate a range of uncertainty of the SRN spectrum arising from the uncertainty associated 
with the model of galaxy evolution. Figure 8 shows the SRN spectra for the SI, S2, II, and 12 models 
for spirals provided that the standard SFR model for ellipticals and a standard cosmological model 
(Slo = 0.2, Ao = 0.8, h = 0.8) are used in common. In order from S2 through SI and then 12 to II, 
supernovae occur more recently (Fig. 4) and the SRN flux becomes higher in the high-energy region of > 10 
MeV, and the SI model gives an intermediate flux. The SRN flux for the II model is about a factor of 2 or 
3 higher than that for the S2 model, and this could be regarded as a maximum uncertainty from evolution 
models of spiral galaxies. On the other hand, the difference among the lower, standard, and higher SFRs 
for ellipticals hardly changes the resulting SRN spectrum. 

Finally, we examine how the SRN spectrum depends on the epoch of galaxy formation, zp. Figure 9 
shows the spectra for zp =3, 4, and 5, using the SI and standard SFR model of galaxy evolution and a 
standard cosmological model (Q,q = 0.2, Ao = 0.8, h = 0.8). Apparently, the dependence on zp is small 
compared to that from the models of spiral galaxies. It is almost negligible above 15 MeV and changes the 
flux at most by a factor of 2 at about 10 MeV. 



4.3. Total SRN Flux 

We calculate the energy-integrated SRN flux in both the whole energy range and observable energy 
range, and compare them with previous estimates. Woosley et al. (1986) estimated the total SRN flux in 
the whole energy range using some different methods, and their estimates spread over a range of 1-250 
cm _2 s _1 . The total SRN flux in our 'standard' model (Slo = 0.2, Ao = 0.8, h = 0.8, zp = 5, the SI model for 
spirals, the standard SFR model for ellipticals, and Cp = 1-93 x 10 8 /i L^QMpc -3 ) gives an moderate value 
of 44 cm _2 s _1 . However, since the SRN cannot be observed in the whole energy range, we have to limit 
the range to avoid various other backgrounds. In the energy range of 15-40 MeV, the other backgrounds 
are quiet (Hirata 1991; Totsuka 1992, and see the next section for more detailed discussion), and we can 
observe the SRN neutrinos through this window, however, only the neutrinos from supernovae which 
occurred in the recent past can contribute to the flux in this range. The SRN flux of our standard model in 
this range is as small as 0.83 cm _2 s _1 , because the supernova rate which is averaged over morphological 
types of galaxies is only less than 0.01 yr _1 in a galaxy with the mass of 10 11 Mq in the recent past, 
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about 10 Gyr after the formation of galaxies (Fig 3). The adopted luminosity density (Efstathiou et al. 
1988) corresponds to the number density of galaxies, ncj = 3.9 x 10 -2 /i 3 Mpc -3 , when all the galaxies are 
assumed to have a mass of IO^Mq, then the supernova rate per unit volume at the present time should 
be less than ~ 3.9 x 10 -4 /i 3 yr -1 Mpc -3 . This value is close to the conservative estimate by Woosley et al. 
(1986) of 2.2 x 10 _4 /i 3 yr- 1 Mpc" 3 , which gives the total SRN flux of only about 1 cm^s" 1 if the rate is 
assumed to be constant throughout the history of galaxy evolution. Most of the total flux of 44 cm _2 s _1 
consists of the neutrinos from supernovae in the early phase of galaxy evolution, which cannot be observed. 
In estimating the SRN flux in Paper I, we used more optimistic values of the supernova rate of 0.1 yr _1 in 
a galaxy and the number density of galaxies ncj = 1.6 x 10 -1 /i 3 Mpc -3 which do not agree with current 
estimates and therefore unlikely. Consequently, the event rate at the SK should also be very small, which is 
the subject of the next section. 



5. Response of the Superkamiokande Detector 

In this section, we calculate the expected event rate at the water Cerenkov neutrino detector, 
Superkamiokande (SK) (Nakamura et al. 1994), and discuss the rate obtained. The procedure of converting 
the SRN flux to an event rate at the SK is described in Paper I. We consider only the reaction v e p — ► e + n 
because the cross section (9.72 x 10 _44 £' e p e cm 2 , where E e and p e are the energy and momentum of the 
recoil positrons measured in MeV, respectively) is much larger than those of the other reactions. The 
detection efficiency of the SK, which has been provided by Totsuka (1994), is taken into account in our 
calculations (see also Paper I). 

We calculate the expected event rate as a function of positron kinetic energy using the SI and 
standard SFR model of galaxy evolution for various cosmological parameters (Fig. 10), and also varying the 
evolution model of spiral galaxies for a standard cosmological model (Slo = 0.2, Ao = 0.8, and h = 0.8) 
(Fig. 11). These figures correspond to Figures 7 and 8, respectively, where the SRN spectra are shown. The 
common feature in all the models is that the event rate is peaked at about 5 MeV with a broad plateau 
over a range of 5-10 MeV, and then decreases rapidly towards higher energies. The event rate does not 
depend strongly on the assumptions of the SFR of elliptical galaxies or the epoch of galaxy formation, zp. 
The event rate of a model with a constant supernova rate is also shown in Figure 10, for which the constant 
rate is normalized such that the total SRN flux is the same as that of the model with Sl = 0.2, A = 0.8, 
and h = 0.8 (thick, middle line). 

Now we consider various backgrounds which may be an obstacle to the SRN detection. They are 
solar and atmospheric neutrinos, electron antineutrinos from nuclear reactors, and non-neutrino terrestrial 
backgrounds (7 rays and [3 decays from unstable nuclei produced by cosmic ray muons or spallation 
products). Below 10 MeV, the SRN will be completely hidden by the reactor i> e 's (Lagage 1985, 1995). In 
the energy range of 10-15 MeV or 10-20 MeV, the main sources of background events are solar neutrinos 
and non-neutrino terrestrial backgrounds, however the background event rate associated with these sources 
drops considerably above 20 MeV(Hirata 1991; Totsuka 1992). Atmospheric neutrinos will dominate the 
SRN above the energy of about 30-40 MeV (Gaisser et al. 1988; Barr et al. 1989). Therefore, we calculate 
the energy-integrated event rate of the SRN in the range of kinetic energy of positrons for 15-40 MeV, 
and also in the more conservative energy range of 20-30 MeV. The energy-integrated event rates in these 
energy ranges for various models of galaxy evolution and the universe, which are considered in this paper, 
are tabulated in Table 2. 
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From the reasons discussed in §4.3, the event rate in the observable range of 15-40 MeV must be quite 
small. When we use the 'standard' model (Slo = 0.2, Ao = 0.8, h = 0.8, zp = 5, the SI model for spirals, 
the standard SFR model for ellipticals, and Cp = 1.93 x 10 8 /i L^Mpc -3 ), the expected event rate at the 
SK is 1.2 yr -1 in such an energy range. The much higher rate of supernova explosions in the early phase 
of elliptical galaxies causes almost no enhancement in the event rate above 15 MeV, since the neutrinos 
emitted from such supernovae are strongly redshifted and fall below the observable energy range. 

The event rate, like the SRN spectrum, shows a rather strong dependence on the evolution model 
of spiral galaxies, the Hubble constant, and the luminosity density of galaxies. It shows only a weak 
dependence on S1 ,A , and the epoch of galaxy formation, zp. Considering all the factors that enhance 
the present supernova rate, we are able to estimate the 'optimistic' event rate. For example, using 
Ho = 100km s~ 1 Mpc~ 1 , zp = 3, the II model for spirals, and the upper end of the B band luminosity 
density C B = (1.93 + 0.8) x 10 8 /i I^Mpc -3 (Efstathiou et al. 1988), the event rate at the SK becomes 4.7 
yr -1 after integrating over the observable energy range of 15-40 MeV. This value is regarded as an upper 
bound of the expected event rate, and if more events are observed, we will have to reconsider a current 
model of evolution of galaxies from which the time-dependent supernova rate is prescribed. 



6. Summary and Conclusions 

We have calculated the SRN spectrum and also the expected event rate at the SK detector by using 
a realistic, time-dependent supernova rate from a reasonable model of galaxy evolution based on the 
population synthesis method. This model of galaxy evolution has been used to investigate the geometry of 
the universe from comparison with the number counts of faint galaxies, and the predicted number counts 
are well consistent with the recent observations. In our analysis the absolute SRN flux in the model is 
determined from the luminosity function and the mass-to- luminosity ratio of nearby galaxies, but not from 
the present supernova rate of our Galaxy which is subject to a large margin of uncertainty. 

In the model we used, the supernova rate is much higher in the early phase of elliptical galaxies, 
and more than half of total supernovae explode during the initial 1 Gyr from the formation of galaxies. 
Thereafter, the supernova rate in spiral galaxies becomes dominant and the total number of supernovae 
until present is consistent with the nucleosynthesis requirement. A constant rate of supernova explosions, 
which was assumed in most of previous papers, is not justified from a viewpoint of galaxy evolution, and 
therefore the SRN spectrum we show here is the most realistic at the present time. 

Supernovae produced at a very high rate in the early phase of ellipticals are responsible for the 
appearance of a hump in the low-energy part ( < 5 MeV) of the SRN spectrum, because the energy of 
neutrinos from those supernovae is degraded by the cosmological redshift effect. This energy degradation 
makes early neutrinos unobservable and the high supernova rate causes almost no enhancement in the event 
rate in the observable energy range of 15-40 MeV at the SK. 

We used a low-density, flat universe with a non-zero A as a standard cosmological model (Sl =0.2, 
Ao = 0.2, and h = 0.8), and investigated the SRN spectrum for other cosmological models relative to the 
standard model. The SRN flux strongly depends on Ho as roughly proportional to Hq. The SRN flux 
changes with Slo and A in a way that smaller Slo and/or larger A suppresses the flux in the high-energy 
region ( > lOMeV), however, these changes are rather small. The shape of the SRN spectrum is almost 
insensitive to either of the above three parameters. 
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We investigated the uncertainty of the SRN spectrum from the model of galaxy evolution. The 
simple, one-zone and n = 1 (SI) model, which we use as a standard evolution model of spiral galaxies, gives 
an intermediate SRN flux among the four different models of spiral galaxies considered in this paper (for 
explanation of the models see §2.1). The difference between the highest flux from the II model and the 
lowest flux from the S2 model is within a factor of 3. In the case of elliptical galaxies, however, the different 
SFRs hardly change the SRN spectrum in either shape or intensity. The effect of changing zp on the SRN 
is almost negligible above 15 MeV, and from zp = 3 to 5 the event rate at the SK changes only by a factor 
of 1.3. 

The total SRN flux in the whole energy range becomes 44 cm _2 s _1 in our 'standard' model (Sl = 0.2, 
Ao = 0.8, h = 0.8, zp = 5, the SI model for spirals, the standard SFR model for ellipticals, and 
C B = 1.93 x 10 8 /i -LbqMpc -3 ), which is a moderate value compared with the previous estimates, however, 
the SRN flux in the observable energy range of 15-40 MeV is as small as 0.83 cm _2 s _1 corresponding to a 
low end of the SRN flux estimated before. This comes from the very small present supernova rate which is 
averaged over morphological types of galaxies with the weight of the mass-to-luminosity ratio. This small 
flux directly leads to a low event rate, and makes the SRN detection at the SK more difficult. 

The event rates at the SK detector are calculated for the various models considered in this paper 
(Table 2), and in the case of the 'standard' model, the expected rate is as small as 1.2 yr _1 in the observable 
range of 15-40 MeV. The event rate depends mainly on the evolution model of spiral galaxies, the Hubble 
constant, and the luminosity density of galaxies, while only weakly on the geometry of the universe (i.e., on 
Slo and Ao) and the epoch of galaxy formation. These properties are in contrast with those of the number 
counts of faint galaxies, and the SRN observation may be complementary to the number counting of faint 
galaxies, if the SRN is detected. From the dependence of the flux on the Hubble constant, the SRN signal 
would provide a new test for measuring the Hubble constant. However, because of presumably low event 
rates, such a test may not be possible in practice. 

Changing the model parameters within the allowable range in a way of increasing the event rate 
(h — ► 1, zp — ► 3, SI— » II, and Cp —> (1.93 + 0.8) x 10 8 /i L B QMpc~ 3 ), we predict the 'most optimistic' event 
rate at the SK, which is 4.7 yr _1 in the energy range of 15-40 MeV. The very small rate of expected SRN 
events at the SK implies the difficulty of gaining some insight on the star formation history in galaxies from 
the SRN observation in the SK experiment. More events exceeding 4.7 yr _1 , if detected, may enforce a 
serious revision in our current understanding of collapse-driven supernovae and evolution of galaxies. 

We would like to thank Y. Totsuka for valuable discussion and useful information about the 
Superkamiokande detector. We would also like to thank P.O. Lagage for the updated data of the reactor 
v e flux. We also appreciate a critical reading of this manuscript by L.M. Krauss. This work has been 
supported in part by the Grant-in- Aid for COE Reaearch (07CE2002) of the Ministry of Education, Science, 
and Culture in Japan. 



- 13 - 



TABLE 1 

Relative Proportions and Mass-to-Luminosity Ratios of Galaxy Types 



Quantity 


E/SO 


Sab 


Sbc 


Scd 


Sdm 


References 


Relative Proportion... 


0.215 


0.185 


0.160 


0.275 


0.165 


1 


{Ml L B ) a 


16.2 


13.0 


9.4 


7.8 


3.4 


2 



a Measured in the B band within the Holmberg radius; in units of h x (M/Lb)q where 
h = Fo/lOOkms-iMpc" 1 . 

References.— (1) Pence 1976; (2) Faber k Gallagher 1979. 

TABLE 2 

Expected Event Rates at the Superkamiokande Detector 



Event Rate (yr r ) 





Ao 


h 


z F 


S-type evolution a 


r i 
Lb 


15-40 MeV c 


20-30 MeV c 


1.0 


0.0 


0.5 


5 


SI 


1.93 


0.63 


0.28 


0.2 


0.0 


0.8 


5 


SI 


1.93 


1.4 


0.62 


0.2 


0.8 


0.5/0.8/1.0 


5 


SI 


1.93 


0.40/1.2/1.9 


0.17/0.51/0.83 


0.2 


0.8 


0.8 


3/4/5 


SI 


1.93 


1.4/1.3/1.2 


0.66/0.54/0.51 


0.2 


0.8 


0.8 


5 


S1/S2/I1/L2 


1.93 


1.2/0.78/1.8/1.5 


0.51/0.34/0.79/0.65 


0.2 


0.8 


1.0 


3 


11 


1.93+0.8 


4.7 


2.1 



a The standard SFR model is used in common to represent the E/SO-type evolution. 

b The luminosity density of galaxies in units of 10 8 /i x Lbq Mpc -3 where h = _ffo/100kms _1 Mpc _1 . 

c Kinetic energy of recoil positrons. 

Note. — We refer the 'standard' model to the case with (£1 , A , h, zp, S-type evolution, Cb)= (0.2, 0.8, 
0.8, 5, SI, 1.93), and the 'most optimistic' model to the case given in the last row with (0.2, 0.8, 1.0, 3, II, 
1.93+0.8). 
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Fig. 1. — The differential number counts of galaxies in the infrared K band in a 1 mag 
interval in a unit area of 1 deg 2 , calculated from the model of galaxy evolution for three 
different cosmological models. The SI model is used to represent evolution of spiral galaxies, 
and the standard SFR model for elliptical galaxies (see text). The symbols represent the 
observed counts taken from various references. The data from the Hawaii observations are 
taken from Gardner, Cowie & Wainscoat (1993), and references of the other observations are 
given in the lower right corner of the figure. The error bars stand for the Poissonian error. 
The selection effects are included in the model curves, using the observational conditions of 
Djorgovski et al.'s survey which are given on the upper left corner of the figure. The thick 
line corresponds to a low density, flat universe with fi = 0.2, A = 0.8, and h = 0.8, the 
dashed line to a low density, open universe with fi = 0.2, A = 0.0, and h = 0.8, and the 
dot-dashed line to an Einstein-de Sitter universe with fi = 1-0, A = 0.0, and h = 0.5. 
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Fig. 2. — The differential number counts of galaxies in the infrared K band in a 1 mag 
interval in a unit area of 1 deg 2 ; the same as Figure 1, but for only a low density, flat 
universe (fi = 0.2, A = 0.8, h = 0.8) with the four different models of evolution of spirals; 
SI (thick line), S2 (short dashed line), II (dot and short-dashed line), and 12 (dot and 
long-dashed line). The standard SFR model for elliptical galaxies is used in common. 
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Fig. 3. — Time variation of supernova rates in a galaxy with the mass of 1O 11 M for various 
galaxy types. The time is elapsed from the epoch of galaxy formation. The dashed lines 
show the rates for the individual galaxy types, and the thick line for the weighted average 
over the types (see text). The SI model is used to represent the evolution of spiral galaxies, 
and the standard SFR model for elliptical galaxies. 
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Fig. 4. — Time variation of supernova rates in a galaxy with the mass of 1O 11 M ; the same 
as Figure 3 but for only the weighted average rates with four different evolution models of 
spiral galaxies; SI (thick line), S2 (short dashed line), II (dot and short-dashed line), and 
12 (dot and long-dashed line). The standard SFR model for elliptical galaxies is used in 
common. 
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Fig. 5. — There is shown the time variation of supernova rates in a galaxy with the mass 
of 1O 11 M ; the same as Figure 3 but for only the weighted average rates for three different 
evolution models of elliptical galaxies; the lower SFR (dot and short-dashed line), standard 
SFR (thick line), and higher SFR (short-dashed line) models. The SI model for spiral 
galaxies is used in common. 
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Fig. 6. — Cumulative numbers of supernovae in a galaxy with the mass of 1O 11 M as a 
function of time elapsed from the epoch of galaxy formation. The dashed lines correspond to 
the cases for the individual galaxy types, and the thick line to the weighted average over the 
types. The SI model is used to represent the evolution of spiral galaxies, and the standard 
SFR model for elliptical galaxies. 
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Fig. 7. — Energy spectra of the supernova relic neutrino background (SRN) for various 
cosmological models. All lines represent the differential number flux of electron antineutrinos 
(z> e 's). The thick lines show the SRN spectrum for a low-density, flat universe with fi = 0.2, 
A = 0.8, and three different values of the Hubble constant; h = 0.5 (bottom line), h = 0.8 
(middle line), and h = 1.0 (top line). The dashed line shows the SRN spectrum for a low- 
density, open universe with fi = 0.2, A = 0.0, and h = 0.8, and the dot-dashed line for an 
Einstein-de Sitter universe with fi = 1-0, A = 0.0, and h = 0.5. The SI model is used to 
represent the evolution of spiral galaxies, and the standard SFR model for elliptical galaxies. 
The epoch of galaxy formation is set to be zp = 5. For an illustrative purpose of comparison, 
the SRN spectrum based on a model with a constant supernova rate is also shown by the 
dotted line, for which the constant supernova rate is normalized in such a way that the total 
(energy-integrated) SRN flux is the same with that for the model with fi = 0.2, A = 0.8, 
and h = 0.8 (thick, middle line). 
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Fig. 8. — Energy spectra of the supernova relic neutrino background (SRN); the same as 
Figure 7 but for only a low-density, flat universe (fi = 0.2, A = 0.8, h = 0.8) with four 
different evolution models of spiral galaxies; Sf (thick line), S2 (short dashed line), If (dot 
and short-dashed line), and 12 (dot and long-dashed line). The standard SFR model for 
elliptical galaxies is used in common, and the epoch of galaxy formation is set to be zp = 5. 
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Fig. 9. — Energy spectra of the supernova relic neutrino background (SRN); the same as 
Figure 7 but for only a low-density, flat universe (fi = 0.2, A = 0.8, h = 0.8) with the 
three different values of the epoch of galaxy formation; zp = 5 (thick line), zp = 4 (dashed 
line), and zp = 3 (dot-dashed line). The Sf model is used to represent the evolution of spiral 
galaxies, and the standard SFR model for elliptical galaxies. 



- 24 - 



l - 




Positron Kinetic Energy [MeV] 



Fig. 10. — Expected event rates at the Superkamiokande (SK) detector for various 
cosmological models, as a function of kinetic energy of positrons produced by the reaction 
Vpjp — > e + n. This figure corresponds to Figure 7 of the SRN spectrum. The thick lines show 
the event rate for a low-density, flat universe with fi = 0.2, A = 0.8, and three different 
values of the Hubble constant; h = 0.5 (bottom line), h = 0.8 (middle line), and h = 1.0 (top 
line). The dashed line shows the event rate for a low-density, open universe with fi = 0.2, 
A = 0.0, and h = 0.8, and the dot-dashed line for an Einstein-de Sitter universe with 
fi = 1.0, A = 0.0, and h = 0.5. The SI model is used to represent the evolution of spiral 
galaxies, and the standard SFR model for elliptical galaxies. The epoch of galaxy formation 
is set to be zp = 5. For an illustrative purpose of comparison, the event rate based on a 
model with a constant supernova rate is also shown by the dotted line, for which the constant 
supernova rate is normalized in such a way that the total SRN flux is the same with that 
for the model with fi = 0.2, A = 0.8, and h = 0.8 (thick, middle line). 
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Fig. 11. — Expected event rates at the Superkamiokande (SK) detector; the same as Figure 
10, but for only a low-density, flat universe (fi = 0.2, A = 0.8, h = 0.8) with four different 
evolution models of spiral galaxies; SI (thick line), S2 (short dashed line), II (dot and short- 
dashed line), and 12 (dot and long-dashed line). This figure corresponds to Figure 8 of the 
SRN spectrum. The standard SFR model for elliptical galaxies is used in common, and the 
epoch of galaxy formation is set to be zp = 5. 



